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Abstract 

Background: Molecular epidemiological studies of HIV-1 in South America have revealed the occurrence of 
subtypes B, F1 and BF1 recombinants. Even so, little information concerning the HIV-1 molecular epidemiology in 
Bolivia is available. In this study we performed phylogenetic analyses from samples collected in Bolivia at two 
different points in time over a 10 year span. We analyzed these samples to estimate the trends in the HIV subtype 
and recombinant forms over time. 

Materials and methods: Fifty one HIV-1 positive samples were collected in Bolivia over two distinct periods (1996 
and 2005). These samples were genetically characterized based on partial pol protease/reverse transcriptase (pr/rt) 
and env regions. Alignment and neighbor-joining (NJ) phylogenetic analyses were established from partial env (n = 
37) and all pol sequences using Mega 4. The remaining 14 env sequences from 1996 were previously characterized 
based on HMA-env (Heteroduplex mobility assay). The Simplot v.3.5.1 program was used to verify intragenic 
recombination, and SplitsTree 4.0 was employed to confirm the phylogenetic relationship of the BF1 recombinant 
samples. 

Results: Phylogenetic analysis of both env and pol regions confirmed the predominance of "pure" subtype B 
(72.5%) samples circulating in Bolivia and revealed a high prevalence of BF1 genotypes (27.5%). Eleven out of 14 
BF1 recombinants displayed a mosaic structure identical or similar to that described for the CRF1 2_BF variant, one 
sample was classified as CRF17_BF, and two others were f]pol/Benv. No "pure" HIV-1 subtype F1 or B" variant of 
subtype B was detected in the present study. Of note, samples characterized as CRF12_BF-related were depicted 
only in 2005. 

Conclusion: HIV-1 genetic diversity in Bolivia is mostly driven by subtype B followed by BF1 recombinant strains 
from the CRF12_BF "family". No significant temporal changes were detected between the mid-1990s and the mid- 
20005 for subtype B (76.2% vs 70.0%) or BF1 recombinant (23.8% vs 30.0%) samples from Bolivia. 
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Background 

Latin America contains the third highest estimated 
number of people living with HIV in the world [1]. The 
HIV epidemic in countries from the Southern Cone is 
characterized by the co-circulation of HIV-1 subtypes B, 
Fl, and BF1 recombinant forms [2-4]. Despite such 
overall similarity, the mosaic patterns of the HIV-1 BF1 
Circulanting Recombinant Forms (CRFs_BF) and 
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Unique Recombinant Forms (URFs BF) display some 
important differences among countries. 

The HIV-1 BF1 epidemics in Argentina [5-9], Chile 
[10], Paraguay [11] and Uruguay [12] are dominated by 
the CRF12_BF and URFs_BF sharing CRF12-related 
structures, in other words sharing some of CRF12 BF 
breakpoints. In addition, new local CRFs_BF with a 
CRF12-related structure have been identified within 
each country, including CRF17 BF in Argentina and 
Paraguay [9,11], CRF38_BF in Uruguay [12], and 
CRF44_BF in Chile [13]. By contrast, the HIV-1 BF1 
Brazilian epidemic is driven mainly by a myriad of dif- 
ferent URFs_BF and at least five distinct CRFs_BF 
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[14-19], unrelated to CRF12BF, which usually show no 
common mosaic recombinant profile among one 
another. Another important difference is that some 
"pure" full-length subtype Fl genomes have been 
described in Brazil [20], while the frequency of such 
full-length subtype Fl viruses or even non-recombinant 
subtype Fl pol sequences is almost absent in other 
South American countries [8,12,21,22]. 

Bolivia is a landlocked country in central South Amer- 
ica, bordered by Brazil to the north and east, Paraguay 
and Argentina to the south, and Chile and Peru to the 
west. The total number of HIV/ AIDS cases reported in 
Bolivia was 5,184 from 1984 through March 2010 [23], 
with an estimated country-wide underreporting rate of 
over 70% [24]. UNAIDS estimates that approximately 
12,000 people in Bolivia were living with HIV in 2009, 
and of those cases, 500 to 1500 were new infections [1]. 
These estimates demonstrate that the Bolivian epidemic 
is still potentially growing. Very scarce data on HIV epi- 
demiology in Bolivia are available, with only two major 
previous studies conducted in the 1990s, which revealed 
the occurrence of subtypes B, Fl and BF1 recombinant 
forms [22,25], similar to what has been described for 
other countries from the Southern Cone [2-13]. In these 
previous studies, however, most samples were geneti- 
cally characterized on the basis of HMA env analysis, 
and very few (n = 14) HIV-1 nucleotide sequences from 
Bolivian seropositive patients are available at the Los 
Alamos HIV database, from those only three were full 
length genomes that reveal the circulation of CRF12 BF 
[5] and CRF17_BF. 

The objective of the present study was to assess the 
trends in the HIV molecular epidemiology in Bolivia 
based on the phylogenetic analysis of HIV-1 pol and env 
sequences at two different time points (1996 and 2005). 

Results 

In order to assess the trends of the HIV-1 molecular 
epidemiology in Bolivia, 62 samples were collected at 
two points in time, 1996 (30 samples) and 2005 (32 
samples), and submitted to amplification and sequencing 
of both the pol (pr/rt) and env genomic regions. Eleven 
samples (nine from 1996 and two from 2005) failed to 
be amplified or sequenced from the proviral DNA in 
one or both genomic regions and were excluded from 
the study. The remaining 51 samples showed a predomi- 
nance of the B subtype (72.5%) and BF1 recombinants 
(27.5%) (Figure 1). Of note, no pure subtype Fl genomes 
(Flpo//Fle«v) were detected. Fourteen highly supported 
(> 95%) monophyletic clusters of at least two sequences 
could be detected in the pol (Figure la) and 11 in env 
region (Figure 2b). Some clusters were composed of 
sequences retrieved from individuals with a known 
direct epidemiological link, such as sexual partners 



(clusters 2a, 5 and 6a), and a mother/child pair (cluster 
7). No apparent direct epidemiological relationships 
were detected for the remaining 10 clusters (1, 2, 3, 4, 6, 
8, 9, 10, 11, 12). Among the 12 monophyletic clusters, 
nine displayed a BpollBenv profile (1-9), two were 
BY\pollY\env (10, 11), and one was Ylpol/Benv (12) 
(Figure 1). 

The 14 BF1 recombinant samples identified in this 
study displayed the following recombination profiles: 
BYlpol/Flenv (n = 11), BYlpol/Benv (n = 1) and YXpoll 
Benv (n = 2) (Figure 1). A more detailed evaluation of 
the 12 BF1 pol recombinant samples was performed 
using bootscan and SplitsTree analyses. From these, five 
samples displayed a recombination profile in pr/rt iden- 
tical to the CRF12_BF, being classified as CRF12_BF-like 
(96BOLP29, 96BOSC14, 05BO06, 05BO15, 05BO18), 
and one sample exhibited a CRF17_BF-like mosaic pat- 
tern (96BOLP23) (Figure 2). 

No epidemiological data were available to confirm the 
proposed close genetic relationship between the 
CRF17_BF (BO02.BO119) strain available in the Los 
Alamos Sequence Database and the BF recombinant 
genome described in the present study. The remaining 
six BF1 pol samples segregate in two highly supported 
monophyletic clusters and exhibited two different 
CRF12_BF-related (Fl/B/Fl/B) mosaic profiles (Figure 
2). In the first cluster (cluster 10), which is composed of 
samples 05BO09 and 05BO20, the first (Fl/B) and sec- 
ond (B/Fl) recombination breakpoints at pol coincide 
with those observed in the CRF12 BF, while the third 
one (Fl/B) is displaced downstream. In the second clus- 
ter (cluster 11), which is composed of samples 05BO03, 
05BO16, 05BO22 and 05BO24, only the third (Fl/B) 
recombination breakpoint at pol coincides with those 
observed in the CRF12 BF, while the first (Fl/B) and 
the second (B/Fl) recombination breakpoints at pol 
were displaced upstream in comparison to those 
detected in CRF12_BF (Figure 3). None of the CRFs_BF 
pol mosaic patterns described in Brazil (CRF28, CRF29, 
CRF39, and CRF40), Uruguay (CRF38), or Chile 
(CRF44) was observed for the HIV-1 samples obtained 
from Bolivia. 

Table 1 summarizes the subtype distribution and 
recombination profiles of the 51 HIV-1 samples from 
Bolivia according to the time period. No significant tem- 
poral change in the HIV-1 subtype distribution was 
observed in Bolivia between 1996 and 2005. During this 
period the prevalence of subtype B decreased from 
76.2% to 70%, while the prevalence of BF1 recombinants 
increased from 23.8% to 30%. 

Discussion 

The present study represents the most comprehensive 
analysis of the HIV-1 molecular diversity in Bolivia 



Guimaraes ef al. Virology Journal 2012, 9:19 
http://www.virologyj.eom/content/9/1 /1 9 



Page 3 of 7 



la 




hah 




96BOSC08 
j 96BOLP22 

1 96BOLP26 ■ 

96BOLP30 
— 96BOSC13 

96BOLP24 . 

11H1 r 96BOSC06 I 

1 96BOSC07 1 

96BO20 1 

96B019 I 




96BOSC03 

0SBO29 
05BO30 
96BOLP2S 

96BOSC1S 
.BR.2002.02BR008 
05BO19 



05BO04 
058005 

96BOSC11 

05BO28 

96BOSC05 

0SB011 

05BO3S h 

05BO37 I 9 

— 05BO14 
05BO23 

-28 BF.BR.1999.BREPM12313 

.20O2.BREPM119 



29 BF.BR.1999.99UFRJ 1 
28 B F.BR. 1999. BREPM 12609 

29 BF.BR.2001.BREPM16704 

39 BF.BR.2004.04BRRJ179 

39BF.BR.2O03.03BRRJ103 




39 BF.BR.2003.03BRRJ327 

17 BF.BO.2002.BO02 B0L119 
96BOLP23 
17 BF.PE.2002.PE02 PCR0155 
- 17 BF.AR.2002.AR02 ARG1139 
— 17 BF.PY.2002.PY02 PSP009S 

17 BF.AR.1999.ARMA038 

38BF1.UY.2004.UY04 4022 

38 BF1.UY.2005.UY05 4752 
38BF1.UY.2O03.UY03 3389 

38 BF1.UY.20O4.UYO4 3987 



-44 BF.CL.2001.CH12 
BF.CL.20qp.CH80 

058009 T 1 n 

05BO20 I 1U 
05BO16 
05BO24 

05B003 



lb 



11 



05BO22 

05B015 

05BO06 

— 12 BF.UY.1999.URTR35 
05BO18 

12 BF.UY.1999.URTR23 

- 96BOSC14 

12 BF.AR.1999.ARMA159 

9GBOLP29 

40BF.BR.2004.04BRRJ115 

40 BF.BR.2005.05BRRJ200 

40BF.BR.2005.05BRRJ055 



40 BF.BR.2C04.04BRSQ46 
r 96BOSC09 " 



96BOSC10 f 

F1.FI.1993.FIN9363 

BR.1993.93BR020 



12 ^ 



CRF28/ 
CRF29 

CRF39 
CRF17 

CRF38 
CRF44 



CRF12 
CRF40 

Fl 

C 



05BO05 



05B031 
B.B0.1996.96BOLP28 

05BO28 

|— 05B029 
058030 



1U0 L 



05BO10 
- 05BO13 



} 1 



} 6a 



0SBO12I 
05BO17 f 



.R.05.BREPM1C81.FJ195091 
D5BOU 



B.BO.1999.BOL0122 

05BO35 

0SBO37 



05BO25 
05B026 

B.AR.99.ARMA132.AY037282 

B.FR.83.HXB2 LAI 1MB BRU.K03455 

B.BO.1996.96BOSC09 

B.BO.1996.96BOSC10 

[ B.B0.1996.96B019 

1 B.6O.1996.96BO20 



}9 

}12 
}5 



05BO14 



.R.01.ARG4006 

— C.BR.98.98BR004 



F1.FR.96.96FR MP411 

Fl.RO.96.BG R07 

F1.AO.06.AO06ANG32 
05BO15 



0SBO20 




J- 10 



F1.B0.1996.96BOLP29 

F1.BR.93.93BR020 1 

F1.BO.1990.96BOLP23 
17 BF.BO.2002.BO02 B0L119 

05BO06 

058018 

058016 

05BO24 
05BO22 
— 05BO03 
-BF.BO.1999.BOL0137 



Fl 



Figure 1 Phylogenetic analysis of 51 HIV-1 pol (a) and 37 env (b) samples from Bolivia Bootstrap values >70% supporting nodes are 
shown on the left of each node and were estimated using the neighbor-joining with Tamura-Nei substitution model. The scale bar indicates 
10% nucleotide sequence divergence. HIV-1 samples analyzed in the present study were in bold format. Arrows indicate those samples with 
intergenic recombination. Highly supported {> 95%) monophyletic clusters are highlighted and numbered. 



performed to date. From a total of 51 samples investi- 
gated in both pol (pr/rt) and env (C2-V3) regions, 72.5% 
of cases were classified as "pure" subtype B, and 27.5% 
were classified as BF1 recombinants. Eleven out of 14 
BF1 recombinants displayed a mosaic structure identical 
or similar to that described for the CRF12 BF strain, 
and the remaining two BF1 recombinants displayed a Fl 
(pol)/B(env) mosaic structure. Subtype Fl {pol)/¥l (env) 
strains were not detected in the present study. 

The high frequency of BF1 recombinants with 
CRF12_BF-like or CRF12_BF-related structures and the 
absence of non-recombinant subtype Fl strains found in 
Bolivia reveal a scenario closely related to that found in 
Argentina, Chile, Paraguay and Uruguay [6-12], and dif- 
ferent from that found in other neighboring countries, 
such as Brazil and Peru [17-19,22]. This pattern suggests 
a more intense HIV-1 influx into Bolivia from the 
southern border (Argentina and/or Paraguay) than from 



the eastern (Peru) or western/northern (Brazil) frontiers. 
It is interesting to note that none of the subtype B sam- 
ples studied presented the typical Brazilian motif 
GWGR in the top of the v3 loop (B" variant of subtype 
B), which is highly (40%) represented in Brazil [26]. The 
reason for such potential unequal viral genetic influxes 
into Bolivia from neighboring countries remains 
unknown and deserves further investigation. 

In addition to the widely disseminated CRF12 BF var- 
iant, other CRFs_BF variants with a CRF12_BF-related 
structure were also described as circulating in Argentina 
and Paraguay (CRF17_BF), Uruguay (CRF38_BF) and 
Chile (CRF44_BF) [9-13]. One sequence with a 
CRF17_BF-like recombinant structure was also detected 
in our study, confirming the circulation of this CRF in 
Bolivia since 1996. Of note, we also detected one highly 
supported cluster (cluster 11) of four BF1 sequences dis- 
playing the same CRF12_BF-related (pol)/¥l (env) 
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Figure 2 SplitsTree analysis of the HIV-1 Bolivian BF1 pol regions. The splits graph was constructed using NeighborNet methodology with 
pairwise distance input, which was estimated by F84 distance. The scale bar indicates 10% nucleotide sequence divergence. Reference 
sequences were indicated by subtype letter or CRF number. The alignment was prepared with HIV-1 group M subtypes B and Fl and CRF_BFs 
from South America. Four clusters containing studied samples were detected and are indicated with a circle. Representative bootscanning plots 
were generated for each cluster. Bootscan analysis was performed on a sliding window of 200 nt of the query sequences moving along an 
alignment of reference sequences by increments of 20 nt. Reference samples used for these analyses are discriminated in the figure legend. 



mosaic profile, obtained from epidemiologically unlinked 
HIV-1 seropositive patients from Bolivia. Such a result 
may represent a new local Bolivian CRF BF variant. 
Full-length genome analysis of these four BF1 sequences 
will be necessary to confirm this hypothesis. 

It is interesting to note that all BF1 recombinants with 
a CRF12_BF-related structure detected in Bolivia dis- 
played larger subtype B segments at pol than CRF12BF. 
A similar finding was also observed in Argentina [5-9] 
and Chile [10]. The generation of CRF12_BF-related 
recombinants with larger subtype Fl fragments in those 
countries is unexpected because of the scarcity (or com- 
plete absence) of pure subtype Fl variants outside of 
Brazil. It has been proposed that the CRF12BF strain 
and other CRFs_BF and URFs_BF variants with related 
structures circulating in South America, as well as other 
groups of related recombinants should be considered 
components of the same "family" of recombinants [3,27]. 



Since its emergence around the early 1980s, the 
CRF12 BF variant has displayed an initial exponential 
expansion in Argentina and Uruguay [28]. Such a fast 
rate of spread and the high prevalence (> 50%) of 
CRF12 BF and CRF12_BF-related mosaic forms in those 
countries could be interpreted as evidence of the high 
viral fitness and/or transmissibility of the CRF12 
"family" of recombinants. However, much lower preva- 
lence (< 25%) of CRF12_BF-like and CRF12_BF-related 
variants has been found in Bolivia [this study], Chile 
[10] and Paraguay [11]. Our results did not depict a sig- 
nificant increase in the proportion of those mosaic HIV- 
1 forms in Bolivia between 1996 and 2005, nevertheless 
all samples from 2005 were collected in only one city: 
La Paz. This finding suggests that the rate of expansion 
and the resulting prevalence of CRF12 BF "family" var- 
iants in South America are country-specific and might 
be dependent on non-virological factors. It has been 
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Figure 3 Schematic representation of the BF1 mosaic structure from clusters 10, 11 and CRF12. Both genomic structures were drawn by 
using the Recombinant Draw Toll available in the Los Alamos homepage (http://www.hiv.lanl.gov/content/hiv-db/DRAW_CRF/recom_mapper. 
html). 



also proposed that the prevalence of CRF12_BF and 
related BF1 strains is higher in countries with more 
extensive intravenous drug epidemics, possibly due to a 
founder effect and the rapid spread of HIV in these 
groups [22]. 

Conclusion 

In summary, our results demonstrate that HIV-1 genetic 
diversity in Bolivia is mostly driven by subtype B 



Table 1 HIV-1 genetic forms identified in the 51 HIV-1 
samples from Bolivia in two distinct time periods (1996 
and 2005) 



Genotype 


1996 (n = 21) 


% 


2005 (n = 30) 


% 


B po//B env 


16 


76.2 


21 


70 


BF1 recombinants 


5 


23.8 


9 


30 


F1 po//B env 


2 


9.5 


0 


0 


CRF1 2-like pol/B env 


1 


4.8 


0 


0 


CRF1 2-like po//F1 env 


1 


4.8 


3 


10 


CRF1 2-related po//F1 env 


0 


0 


6 


20 


CRF1 7-like pol/FI env 


1 


4.8 


0 


0 



List of abbreviations: HIV-1 Human immunodeficiency virus-1, pol Polymerase, 
pr/rt Protease/reverse transcriptase, env Envelope, NJ neighbor-joining, CRF 
Circulating recombinant forms, URF Unique recombinant forms, AIDS Acquired 
immune deficiency syndrome, HMA Hetereduplex mobility assay, DNA 
Deoxyribonucleic acid 



followed by BF1 recombinant strains from the 
CRF12BF "family" and no significant temporal changes 
in the subtype distributions were detected between the 
mid-1990s and the mid-2000s. 

Materials and methods 

Study population 

A set of 62 samples were collected in Bolivia in 1996 (30 
samples) and 2005 (32 samples). Demographic, labora- 
tory and clinical data of Bolivian subjects collected in 
1996 were previously described [24]. Samples from 1996 
were collected in three distinct Bolivia cities [La Paz (n 
= 9), Cochabamba (n = 5) and Santa Cruz (n = 16)]. 
The second cohort was composed exclusively of patients 
living in La Paz, and most of them were under antiretro- 
viral therapy. All participants signed the informed con- 
sent form, and the study was approved by the Bolivian 
Ministry of Health. 

HIV-1 DNA amplification and sequencing 

Genomic DNA was extracted from whole blood by the 
phenol/chloroform method from samples collected in 
1996, as described [25]. For samples collected in 2005, 
the QIAamp DNA kit was used (Qiagen Inc., CA, U.S. 
A.) according to the manufacturer's protocol. 
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Amplification of HIV-1 env and pol regions was per- 
formed using an in-house nested PCR protocol method 
as described elsewhere [29]. Samples collected in 1996 
were previously subtyped based on the HMA env-gpl20 
C2-V3 region and 11 of them were also sequenced and 
phylogenetically analyzed as described [25]. Partial pro- 
tease/reverse transcriptase (pr/rt) sequences were gener- 
ated for the samples from both periods, and new env 
sequences were obtained only from samples collected in 
2005. No HIV-1 seropositive individuals from the first 
time period were undergoing antiretroviral therapy at 
the time of blood collection; however, more than a half 
of the individuals in the second group were being 
treated. 

HIV-1 subtype classification 

HIV-1 subtype classification was performed by phyloge- 
netic analysis with the neighbor-joining method using 
the Tamura-Nei model as implemented in Mega 4.0 
[30]. Those samples that present divergent clustering in 
the phylogenetic analysis were investigated for recombi- 
nation by bootscan analysis as implemented in the Sim- 
plot version 3.5.1 (http://www.med.jhu.edu/deptmed/ 
sray/download/), using representative HIV-1 sequences 
from group M. Bootstrap values supporting branching 
with reference sequences were determined in neighbor- 
joining trees constructed using the Kimura 2-parameter 
model, based on 100 re-samplings, with a 200 nucleotide 
sliding window moving in steps of 20 bases. The Split- 
sTree program version 4 [31] was employed to confirm 
the phylogenetic relationship of the recombinant sam- 
ples using the NeighborNet based on the pairwise dis- 
tance estimated by F84 parameter model. 

Sequence data 

GenBank accession numbers: JN710748-JN710798 (env- 
gpl20 C2-V3 region), JN710799-JN710830 (pr/rt). 
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